hep-th/0502227 

CERN-PH-TH/04-251 

IFT/05/11 

MIT-CTP-3591 

TIFR/TH/05-03 



in ■ 
o 

Finite temperature effective action, A0IS5 black holes, 

and 1/N expansion 

00 ■ 
(N ; 

Luis Alvarez- Gaume^, Cesar Gomez 11 , Hong Liu*, and Spenta R. Wadia 1 ''* 

r-> ■ 
<n : 

CN ■ Physics Dept. -Theory Division, CERN, CH- 1211, Geneva 23, Switzerland 

<N ■ 

^Institute de Fsica Terica, Universidad Autnoma de Madrid, Cantoblanco 28049 Madrid, Spain 

in 

* Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 



'Dept. of Theoretical Physics, Tata Institute of Fundamental Research, Mumbai, 400 005, IndiM 



We propose a phenomenological matrix model to study string theory in AdS 5 x S5 in 
the canonical ensemble. The model reproduces all the known qualitative features of the 
theory. In particular, it gives a simple effective potential description of Euclidean black 
hole nucleation and the tunnelling between thermal AdS and the big black hole. It also 
has some interesting predictions. We find that there exists a critical temperature at which 
the Euclidean small black hole undergoes a Gross- Witten phase transition. We identify 
the phase transition with the Horowitz-Polchinski point where the black hole horizon size 
becomes comparable to the string scale. The appearance of the Hagedorn divergence 
of thermal AdS is due to the merger of saddle points corresponding to the Euclidean 
small black hole and thermal AdS. The merger can be described in terms of a cusp (As) 
catastrophe and divergences at the perturbative string level are smoothed out at finite 
string coupling using standard techniques of catastrophe theory. 
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1. Introduction 



The AdS/CFT correspondence has enabled us to begin understanding various aspects 
of quantum gravity in a more quantitative way [[[]. In particular reliable computations of 
black hole thermodynamics have been possible using the AdS^/CFT correspondence 
In this paper we would like to address some aspects of AdS 5 black hole physics in the 
context of type IIB string theory in AdS 5 x S 5 , using the dual gauge theory at finite 
temperature. 

The thermodynamic aspects of quantum gravity in AdS spacetime were discussed long 
ago in an important paper by Hawking and Page |J, who realized that it is possible to 
define a canonical ensemble for quantum gravity and in particular for Schwarzschild black 
holes. They found that an asymptotic AdS spacetime allows two Schwarzschild black hole 
solutions, which are since called small black hole (SBH) and big black hole (BBH). As the 
names suggest SBH can have a horizon radius that can be very small compared to the 
size of AdS and BBH has a horizon radius which is comparable to (or much larger than) 
the radius of AdS$. Further SBH has negative specific heat and is unstable, while BBH 
has positive specific heat and is stable (meta-stable.) Hawking and Page also found that 
the system undergoes a first order phase transition at a temperature T\ comparable to the 
inverse curvature radius of the spacetime. Below Ti, the system is described by a thermal 
gas in AdS, while above T\ it is described by a BBH. With the discovery of the AdS/CFT 
correspondence Witten realized that a BBH in AdS 5 is naturally described by 

the deconfinement phase of M = 4 Super- Yang-Mills (SYM) theory on^xS 1 . He argued 
that the Hawking-Page transition corresponds to a large N deconfinement transition in 
the gauge theory at strong coupling. 

Several authors P,P,|TUf have discussed the partition function of the free M — 4 SYM 
theoryi and found that the large iV deconfinement transition persists at zero coupling. In 
particular it was found that the deconfinement transition happens exactly at the Hagedorn 
temperature of the low temperature thermal AdS phasei. Near the Hagedorn temperature, 
the free energies of both high and low temperature phases become divergent and string 



2 See [11,1^,13,14,15,16,17,18] for other recent discussions of phase transitions in weakly cou- 
pled Yang-Mills theory. 

3 At strong coupling, the Hagedorn temperature for the thermal AdS is much higher than the 
Hawking-Page temperature. 
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perturbation theory! breaks down. In [19| the smoothening of the Hagedorn transition 
at finite string coupling was discussed. This requires a careful understanding of nonper- 
turbative effects in 1/N. It was found that the divergences in perturbation theory are 
removed by two distinct mechanisms. The divergent terms in the high temperature phase 
can be resummed, leading to a noncritical string description. This happens at a scale 

4 

T — Th ~ N~3 . The Hagedorn divergence of the low temperature phase is removed by 
summing over the contributions from the thermal AdS and the noncritical string back- 
ground, happening at T — Th ~ N~ 2 . 

In this paper we extend the analysis of fn| to finite £ t Hooft coupling and study 
various non-perturbative aspects of black hole physics in AdS 5 using the boundary gauge 
theory. The latter is precisely formulated but technically difficult to deal with in the strong 
coupling region, where one can make contact with gravity. In such circumstances one is, 
naturally led to an effective action approach, relying on the choice of an order parameter 
and its symmetries. The difficulty is then transferred to the coupling and temperature 
dependence of coefficients of the effective action coding the microscopic theory. In spite of 
this difficulty one can hope to make progress. Above all, one is encouraged by the success 
of a similar programme in QCD. 

The strategy has two parts. First one may try to extract certain universal features of 
string theory in AdS 5 from the effective action. The hope is that universal features do not 
depend on the exact details of the effective action and can be extracted by exact analysis 
of a tractable model. Secondly one can approximately determine the coefficients of the 
effective action by explicitly matching, in our case, with data in the dual supergravity 
description. 

String theory backgrounds like thermal AdS, BBH and SBH appear as saddle points 
in the Euclidean path integral of Yang-Mills theory. Perturbative string expansion around 
each of them is given by the large iV expansion around the corresponding saddle point 
in the boundary theory!. As one varies the temperature, such expansions break down at 
various places where their coefficients develop nonanalytic behavior. One example is the 
Hagedorn temperature of thermal AdS. The other is the temperature (called T by [|J) at 



4 Since the 1/N expansion in the boundary theory corresponds to the perturbative string 
expansion in the bulk, in this paper we will use the word "large iV expansion" and "perturbative 
string expansion" interchangeably. 

5 For the perturbative 1/N expansion around the saddle point, it is not important whether 
the saddle of interest is stable, metastable or unstable. 
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which SBH and BBH saddles appear. While such nonanalytic behavior is rather puzzling 
and hard to understand from the perturbative string point of view, in the dual Yang- 
Mills theory they arise due to the fact that in the large N limit the number of degrees 
of freedom goes to infinity. The non-analyticity occurs for the same reason as in the 
thermodynamic limit of classical statistical physics. At these temperatures, as in the case 
analyzed in |19[ , a non-perturbative treatment is required no matter how large N is. Thus 
these non-analyticities are excellent probes of the non-perturbative structure of the theory. 
Moreover, since their appearance is intrinsically tied to the large iV limit, it is expected 
thnat they possess a certain degree of universality, just as in critical behavior in condensed 
matter systems. One may hope that the qualitative behavior at these critical temperatures 
sould be insensitive to the precise details of the theory and could be captured by studying 
much simpler systems. This gives us hope that we can study critical behavior, and hence 
non-perturbative aspects of large iV Yang-Mills theory at finite or strong coupling, and in 
turn yielding insights into string theory in AdSs. 



It was discussed in [[10] that the partition function of the SYM theory can be written 
as a matrix integral over the effective action of the Polyakov loop. For Yang-Mills theory at 
finite coupling, we have no way of computing this effective action explicitly. Nevertheless, 
with universality in mind, here we propose a class of effective actions as "phenomenological 
models" to approximate the full theory. We show that models in the class have a large N 
phase structure resembling that of a weakly coupled string theory in AdSs x S$. This gives 
strong indication that strongly coupled Af = 4 SYM theory belongs to the same universality 
class. This also gives us reason to believe that critical behaviors of the bulk string theory 
at places where string perturbation theory breaks down can indeed be captured by much 
simpler models. 

The simplest model within the class, which we will refer to as (a, b) model, can be 
considered as a truncation of the full effective action of the theory to the lowest nonlinear 
termsi. Being exactly solvable to all order in N, this model provides an ideal representative 
to study the critical behaviors of the universality class. We proceed to perform a detailed 
study of various non-perturbative aspects of this model. The results, when translated into 
the language of bulk string theory, can be summarized as follows: 



6 The model contains two parameters (a, 6), both are functions of the 't Hooft coupling A 
and temperature T. We will assume some qualitative dependence of (a, b) on T as part of the 
phenomenological input data. This model has been briefly discussed earlier in the mean field 



approximation as a toy model for weakly coupled gauge theories in [10]. 
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1. We give an effective potential description of the tunnelling between thermal AdS and 
the BBH. The Euclidean SBH plays the role of the bounce (also called a thermalon). 
We compute the tunnelling rate in our effective theory. 

2. We find that the Euclidean SBH undergoes a third order Gross- Witten pO| , pT| , p2| 



phase transition in the large iV limit at a temperature T c below the Hagedorn temper- 
ature. We identify the phase transition with the Horowitz-Polchinski correspondence 



point []23[ where the event horizon size of the SBH becomes comparable to string scale. 

3. The breakdown of the perturbative string expansion of thermal AdS at the Hagedorn 
temperature is due to the merger of the saddle points corresponding to the SBH and 
thermal AdS. The merger can be described in terms of a cusp (A3) catastrophe. The 
simplest possibility allowed by the symmetry. Similarly the breakdown of perturbative 
string theory around the BBH when it merges with the SBH saddle can be understood 
in terms of a fold (A2) catastrophe. The divergences at the perturbative level are 
smoothed out at finite N using the standard techniques of catastrophe theory. 

4. A common theme in our study of the critical behavior when a CFT approaches a 
singular point is that there always exists a double scaling limit and it is likely that 
the theory in the double scaling limit is described by a noncritical string background. 



This also resonates with the result of [19| and the behavior of other singular CFTs 



discussed in 24,25,26 



While these features are studied explicitly only in the simplest (a, b) model, we believe 
they persist for all models in the class due to universality of the large N phase transition 
and the catastrophe. 

The plan of the paper is as follows. In the next section we review some aspects of the 
thermodynamics of quantum gravity in AdS*, x £5 which we aim to reproduce in the large 
N limit of our "phenomenological" models. In section 3, we review some aspects of the 
computation of the Yang-Mills partition function using the effective action of the Polyakov 
loop and present the truncated models. Section 4 is devoted to a detailed study of the phase 
structure of the (a, b) model at large N. We discuss in detail the thermal history of the 
theory in the canonical ensemble. We also show that the sharp Hawking-Page transition in 
supergravity is smoothed out to a finite cross region at finite N. In section 5 we elucidate 
the role of SBH as the bounce which mediates the tunnelling between BBH and thermal 
AdS (and vice versa depending on the temperature) and calculate the tunnelling rate. We 
also connect the bounce and the large order behavior of perturbative theory. In section 
6 we study the critical behaviors of the theory at temperatures where the perturbative 
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string expansion breaks down around at least one of the three backgrounds. They can 
be understood using catastrophe theory. We conclude in section 7 with a discussion of 
future directions. We have also include a few appendices which contain details of some 
calculations. 

2. Review of Hawking-Page transition in Euclidean Quantum gravity 

In this section we review the results of , to be reproduced using the matrix models 
in later sections. 

The canonical ensemble for quantum gravity in AdS can be defined as a path integral 
over the metric and all other fields asymptotic to AdS with time direction periodically 
identified with a period [3 = 1/T. At semi-classical level, i.e. R 2 /lp ^> 1, where R is 
the curvature radius of AdS, such a path integral is dominated by configurations near the 
saddle points, i.e. classical solutions to the Einstein equations. If we assume spherical 
symmetry and zero charge, there are three possible critical points, which are thermal 
AdSs (Euclidean AdS with time direction periodically identified), a big (Schwarzschild) 
black hole (BBH) and a small black hole (SBH). Among them thermal AdS and BBH 
are locally stable, while SBH has a negative mode and it is unstable. The thermal AdS 
background has topology S 1 x i? 4 , while SBH and BBH have topology R 2 x S 3 , all of them 
with a common boundary S 1 x S 3 . The Euclidean time direction in black hole backgrounds 
are contractible and the winding numbers are not conserved. In contrast the time circle 
in thermal AdS is noncontractible and the winding number is conserved. 

The classical action for thermal AdS is I\ = 0. This is standard in string theory: with 
a noncontractible time circle, there is no genus zero contribution to the free energy. A 
Schwarzschild black hole solution exists in AdS only for a Hawking temperature greater 
than _ 





(2.2) 



The corresponding classical Euclidean action is given by 




(2.3) 
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where 2k 2 is the five-dimensional Newton's constant!. We will denote the classical 

actions for large and small black hole respectively. The specific heat of the large black 
hole is positive and thus it is thermodynamically stable (i.e. it can reach locally stable 
thermal equilibrium with thermal radiation). The small black hole has a negative specific 
heat. The action 7_ of the small black hole is always greater than the action of thermal 
AdS and of the big black hole. At temperature 

r > = ia > T ° (2 ' 4) 

the action for the big black hole is 1+ = = I±. When To < T < T\, i+ > 0, and the 
saddle corresponding to thermal AdS dominates. When T > T\, 1+ < 0, the big black 
hole (BBH) dominates. There is a change of dominance at T\. This is the Hawking-Page 
transition. In the classical limit k 2 — > 0, this is a sharp first order transition. We expect 
that at finite k 2 the transition should be smoothed out. This we will see explicitly in the 
gauge theory description. 

In the Minkowski description the spectrum of fluctuations around AdS 5 and BBH (all 
of positive frequency in Euclidean and Minkowski descriptions) can be interpreted in terms 
of physical particles and they constitute the meta-stable thermal ensemble around these 
backgrounds. The SBH on the other hand has a negative eigenvalue in the spectrum of 
small fluctuations in the Euclidean description. Hence the SBH fits the description of an 
instanton relevant for the tunnelling between thermal AdS$ and the BBH. For example, 
one expects the rate for a BBH to tunnel into a thermal AdS is expected to be 

T x = A ie - {I -- I+) (2.5) 

That is, through thermal fluctuation, a big black hole can turn into a small black hole. The 
small black hole (since it has a negative specific heat) then can either shrink to thermal AdS 
by emitting thermal radiation or grow back into a big black hole by absorbing radiation. 
Similarly, the thermal AdS background has also a nonzero probability to nucleate a small 
black hole with probability 

T 2 = A 2 e~ I ~ (2.6) 
Afterwards the small black hole can shrink back to the thermal AdS or grow into a big 



black hole. The prefactors A%, A 2 in (|2.5| ) and (|2.6| ) are given by the determinants of small 



7 Note that °^ N 2 , where ./V is the rank of the gauge group in the boundary Yang-Mills 
theory. 
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fluctuations around the relevant background. In thermal equilibrium, the probability to 
go from a typical state in thermal AdS to that of a big black hole or back should be the 
same. 

When Tq < T < Ti, the big black hole phase is metastable, since it has a higher free 
energy than that of that of thermal AdS and Fi > 1?2. But string perturbation theory 
around it is well defined until To is reached where we expect the perturbation theory to 
break down. Similarly, when T > T±, thermal AdS becomes metastable and 1?2 > IV For 
a large AdS with R ^> l s (Z s is the string length) the perturbation theory around thermal 
AdS breaks down at a much higher Hagedorn temperature Tjj ~ A. In the Hawking-Page 
discussion, there also exists a temperature beyond which the thermal graviton gas in 
AdS will collapse into a big black hole. For a weakly coupled string theory in AdS*, xS$, T2 
is of order — ^-j- and is much higher than the Hagedorn temperature Tjj ~ y- for thermal 
AdS. 



3. Effective action at finite temperature 

In this section we will introduce a phenomenological matrix model for understanding 
string theory in AdSs x £5 at finite temperature. 

We first give some general discussion of the partition function of = 4 SYM theory 
on S 3 . We consider the theory in the canonical ensemble, i.e. the Euclidean time direction 
is periodically identified with a period of f3 = h. It was pointed out in |I(J (see also pTfl ) 
that the Yang-Mills theory partition function on S 3 at a temperature T can be reduced to 
an integral over a unitary U(N) matrix U, which is the zero mode of Polyakov loop on S 3 , 

Z{\T) = J dUe s *ff^ (3.1) 

with 




U = Pexp i / Adr (3.2) 



where A{j) is the zero mode of the time component of the gauge field in S 3 . This follows 
from the fact that apart from A all modes of M = 4 SYM on S 3 are massive. Hence one 
can integrate them out to obtain an effective action for A. Gauge invariance requires that 
the effective action must be expressed in terms of products of trt/ n with n an integer, since 



these are the only gauge invariant quantity that can be constructed from A alone. S e ff(U) 
has a Zn symmetry 

due to global gauge transformations which are periodic in the Euclidean time direction up 
to Zn factors. A generic term in S e ff(U) will have the form 

trU ni trU n2 ■■■trU n \ n x + ■ ■ ■ n k = (mod N), k > 1 

We can expand S e f / in terms of a complete set of such operators, with the first few terms 

S eff (U) = a^vUtrU- 1 + b ± (tvUtvU^) 2 + a 2 tvU 2 trU- 2 + atrUHrUhitf + ■■■ (3.3) 

The coefficients in the expansion are functions of 't Hooft's coupling A, and T. While these 
coefficients are in principle calculable at weak coupling, the explicit computations are in 
general very complicated (see e.g. [p8| ). At finite or large 't Hooft coupling, there is no 



available tool at the moment to attempt such a computation. Even one were able to find 
the expansion ( |3.3|) explicitly, to perform a finite N computation of the matrix integral 
( p.l|) is still a daunting task, if not impossible. 

In order to make progress, in this paper, we will consider the truncation of ( |3.3|) to 
terms containing only powers of tr£7tr£7 -1 , i.e. we consider an effective action of the form 

S eff (U) = S(x), x=^-tvUtvU^ 

b 2 C 3 (3 ' 4) 

= atrUtrll' 1 + (trUtrU^) + j-^ (trl/trU^) H 

Our consideration is phenomenological, motivated by the AdS/CFT correspondence to 
search for effective actions which lie within the same universality class as that of the SYM 
theory at finite coupling. At a heuristic level, one may also consider (|3.4| ) as arising from 
( p.3| ) by "integrating out" all the higher moments trC7 n , trU~ n (n> 1). In [27| it was also 



argued based on group theory considerations that terms of the form (|3.4|) dominate over 
other types terms in ( |3.3|) in the large iV limit. 

We will consider a class of matrix models of the form ( |3.4j ) satisfying the conditions 
that S(x) is convex and S'(x) is concave. We show in the next section and in Appendix 
A that in this class the large iV phase structure appears to be universal. In particular, 



the phase structure precisely reproduces^ the phase structure of a weakly coupled string 
theory in AdS. We believe this is a strong indication that strongly coupled TV = 4 SYM 
theory also lies in this class. The simplest model in this universality class is given by the 



first two terms in (|3.4| ) 



Z(a,b) 



dU exp 



atxUtxU ] + —2 (trUtrU^) 



(3.5) 



with b > 0. 

As discussed in the introduction, we are interested in the critical behaviors in regions 
of parameter space where large N expansions around various saddle points break down. 
Being exactly solvable at finite N, (|3j|) provides a simple, nice representative to study the 
critical behaviors of the universality class. 



To conclude this section, we note that (|3.5| ) was discussed in [|10| in the Hartree-Fock 
approximation. It was noted that for b > a first-order phase transition resembling the 
Hawking-Page transition occurs. This observation was a motivation for the investigations 
in this paper. 



4. Large N phase structure of the universality class 



In this section we study the large iV phase structure of (|3.4| ) and (|3.5| ) (b > 0). For 



our later purpose of studying the critical behaviors of Q3.5| ), we give a detailed discussion 
of the phase structure of (|3.5|) using a method suitable for finite iV analysis. We point out 



general matrix model ( |3.4| ) 



Z=jdUe N2s{x \ x = -^trUtrtft (4.1) 



has the same large N phase structure as (|3.5|) provided that S(x) is convex and S'(x) is 
concave. For completeness we have included in Appendix A an alternative discussion of 
the phase structure of (|4.1|) using the Hartree-Fock method. 



8 We will assume some qualitative dependence of S(x) on T as part of the phenomenological 
input data. 
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4-1- Effective potential 

For b > 0, equation (|3.5|) can be rewritten using a Lagrange multiplier /i: 

TV f°° 2 2 f 

Z(a,b) = / d^ie-^^- a)2 \ dU exp [/jtrLW*] . (4.2) 



The matrix integral in (|4.2| ) can be further simplified by introducing another Lagrange 
multiplier g. For example for //0, one finds 



exp 

(4-3) 



2^ Jo 
with 

..o . r r Mn . 

(4.4) 



e N2F ^=JdUexp 



^(trd + trC/t) 



The formula for ^ < is obtained by taking g — > zg. 

The large iV expansion of (|4.4| ) and the corresponding third order phase transition is 
well known |2g,gl],|29| 



2 

^j- + nonperturbative (7 < 1 or g imaginary 

F(g) = { ' (4.5) 

g-l\ogg-^ + 0(l/N 2 ) g>l 



The order parameter of ( |4.4[ ) can be taken to be 

1 1 r) F 

Mp) = ^<w) s = -<W> s = ^ 



+ ■ ■ ■ g < 1 or g imaginary (4-6) 



characterizing the eigenvalue distribution of U. When < p\ < \ (g < 1), the system is in 
a phase whose eigenvalue distribution does not have a gap on the unit circle. In particular, 
for (7 = 0, the distribution is uniform. When for 1 > pi > | (g > 1), the distribution 
develops a gap. ((L5|) and (|4.6|) do not apply to g ~ 1, where the system undergoes a third 
order phase transition in the large iV limit. At finite iV the third order discontinuity 
in ( |4.5| ) is smoothened out by non-perturbative effects. They will be discussed in later 
sections when needed. 
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One can now rewrite (3J3) as a two dimensional integral 



Z 



N 3 



4V7T& 



dfx 



— / gdge 



-N 2 V(ti,g) 



(4.7) 



with 



V{n,g) 



2 

4 /, 

2 

4 n 

{±- b ^-a) 2 + £-g + \\ogg+l + 0{l/N*) (i > 0, g > 1 



H < 

At > 0, < g < 1 



It is often convenient to integrate out g to reduce ( |4.7|) to a one dimensional integral 

N 



Z(a,b) 



with 



2Vivb 



1 



dfj, e 



-7V 2 Q( M ) 



Q(p) = ^(a^-«) 2 -f(p) 



(4.9) 



(4.10) 



and J-(jj) was defined in (|4.3|) . The large expansion for J-(fj) was found in |I9|] , e.g. the 
leading order terms are 



0-^ log(l - n) + 



/U < 1 



1 w 



21 _ w + ±log(l-w) + 0(^) m> 1 
where for /u > 1 we have introduced 



(4.11) 



IV 



1 



1 



(4.12) 



Inherited from (O), (|4.8| ) has a third order discontinuity at p = 1 which needs to 
be supplemented with a non-perturbative treatment. (|4.10| ) with JF given by ( |4.11| ) has 
divergences^ and first order discontinuity at \x = 1. Again a non-perturbative treatment 
is necessary, as discussed in detail in []19|. Part of the subtlety at fx = 1 in (|4. 10| ) has to 
do with that at li = 1, g becomes massless and the effective potential V(n,g) is flat in 
the range < g < 1. Thus near (j, = 1 it is more convenient to use the two dimensional 
effective potential ( |4.8|) . Also in ( 4.11 ) we have suppressed a subdominant term which 
should be taken into account in the analysis of the phase structure. This is automatically 



The second and higher derivatives of J-(li) are also divergent for li — >• 1+. 
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taken care of in the two-dimensional integral (|4.7|) . In this and following sections we will 
use both forms of the effective potential (|4.8|) and ( f4.10|) depending on convenience. The 
one-dimensional effective potential Q(n) is easier to visualize than the two-dimensional 
potential V(/j,,g). But as we mentioned, equation ( fOf ) is more convenient around fx = 1. 

We note that the trick used in (|4.2| ) to reduce ( |3~5| ) to an integral transform of ( |4~3| ) 
can be generalized to find the large N phase structure of general matrix models (|4.1| ) with 
S(x) convex. Since S(x) is convex, it admits a Legendre transform^: 

S(fx) = Max x (fxx — S(x)) 

The Legendre transform is involutive. If we do it twice we get back the same function. 
Hence S(fi) is also convex. We can then write ( |4.1| ) as 

Z = J dU e N2s ^ =JduJdLt e N2 ^~ s ^ (4.13) 

and the second integral over \x is carried out using saddle points. For large N this will give 
an excellent approximation. Convexity in fact guarantees that there is a unique saddle 
point contributing. If we now exchange the order of integration in ( |4.13| ) and use QL3| ) , we 
find that 

Z = J dLie- N2Q ^ (4.14) 

where 

q(/0 = s(aO—f(/0 ( 4 - 15 ) 

and is given by ( [4 . 1 1| ) in the large iV limit. We will show in next subsection that 

Q(li) leads to the same large iV phase structure as ( |4 . 1 U| ) provided S'(n) is also convexS. 

4-2. Phase structure 

In the large iV limit the critical points of V in ( [4.71 ) describe different phases of the 
theory which in turn correspond to different bulk string theory geometries. The minima 
correspond to (meta)stable phases, while saddle points^! (or maxima) to unstable phases. 
Note that in the large N limit, the eigenvalue distribution of the Polyakov loop U at a 
critical point follows from (|4.4[) with g given by its value at the critical point. Since there 



10 For general properties of Legendre transformations see for instance the book [ ]30| |. 

11 That S'(fi) is convex means that S'(x) is concave. 

12 By saddle here we refer to saddle points of V(li, g) on the real li — g plane. 
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is a one-to-one map ( }4.6|) between p\ and g, g can be considered as the order parameter of 
the theory. After integrating out g one can also interpret \x in ( [4.9| ) as the order parameter 
(at least in the range /xl). 

Before discussing the critical points of the theory in detail (which is somewhat in- 
volved), we note that qualitative features of the critical point structure of ( |3.5| ) can be 
conveniently visualized by plotting the one dimensional effective potential ( |4.iq ) in 

the large iV limit. Depending on the values of (a, b), Q can have one or three critical points 
(see fig. 2). The critical point structure in the (a, b) plane is plotted in fig. 1 (see below). 
Below curve I in fig. 1, Q has one minimum. Between curve I and curve H, it has three 
critical points \i\ < \ii < ^3, with two minima (fii and ^3) and one maximum (1^2) ■ The 
two minima change dominance on curve II. On curve I, \i<i and fis merge. At curve H, (jli 
and (12 merge. To the right of curve H, in addition to ^3, (|4.8| ) also has a tachyonic saddle 
which is not visible in the leading order Q-plot here. 




-1 . 5 



-1 



-0.5 0.5 1 

Fig. 1: This figure plots the critical point structure of the theory in the a — b 
plane. Below line I, there is one critical point. There are three critical points 
between line I and line H, two minima, one maximum. At line II, two minima 
exchange dominance. 



We now describe the critical points of (|4.8|) in detail: 
1. From the first two lines of (|4.8|) one finds the following critical point 



^1 



91 = 



(4.16) 
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Fig. 2: Plot 2a corresponds to (a, b) below curve I of fig. f, 2b to (a, b) between 
curve I and curve II, 2c to (a, b) between curve II and curve H, and 2d to (a, b) 
lying to the right of curve H. Note the discontinuity in the first derivative at fj, = 1 
is due to the large iV approximation. It is smoothened out by non-perturbative 
effects. 



with 

V{n ugi ) = (4.17) 

and 

v "=(l: v vz)=\{\ 4) (418) 
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where V gg denotes and so on. For a < 1, it is a local minimum. V" becomes 
singular for a = 1 and tachyonic for a > 1. Since gi = 0, it describes a uniform 
eigenvalue distribution in the unit circle. 
2. For a < 1 and c = 2 ^ 1 ~°^ < i ? there is an additional saddle point from the second line 
of equation (f4.8|) at 



/* 2 = 1, S2 = v^, c=^ T ^<l (4.19) 

with 

(1 - a) 2 

Vbi,9i) = L 1 £ L (4-20) 



and 

T/" = , 

2 V ~Vc 

14 



y// = 1 ( I +c V°\ (4.21) 



Note that V" has a negative eigenvalue. Since pi(n2i 92) = < h, P2 describes a 
gap less phase in the eigenvalue distribution of U. As a — > 1, ( |4.16| ) and ( |4 . 1 9| ) merge, 
after which ( |4.19| ) disappears while ( [4.16| ) becomes tachyonic. 
3. From the third line of Ql8|) the equations for the critical points are given by 



1 

2b 



{H - a) 



g 



0, 



9_ 
2/x 



1 + 



1 







(4.22) 



4p 2 ' 2p ' 2g 

Note that the second equation in ( f4.22|) only has real solutions for \i > 1, in which 
case one finds that§ 

9 = 1 ~TTT' (4-23) 



1 — w 

where w was introduced in ( [4.12|) . The eigenvalue distribution for such a g is given 
by 

* = 1 * =i + 2£>l (4.24) 



Substituting (^4.23| ) into the first equation of ( |4.22|) one finds an equation for /x 



fx — a 



(w + 1) 



2 b 4 ( 4 ' 25) 

Of course (|4.25|) can also be obtained directly by extremizing the leading order term 
of the effective potential (|4.10|) for (j, > 1, namely 



Q(p) = V(ji,g(ji)) 
Note that 



1 w 



46 



21-w 



2 1 

■7 + ^7 + 



Ilog(l-«,)+o(^ 



b {I 2 y/JT^l^/ 2 



(4.26) 



(4.27) 



Deriving \i from ( |4.25j ), g and pi can then be found from ( 4.23 ) and ( 4.24 ). Since pi 
is a monotonic function of fj,, one can treat Q((J.) as an effective potential for p\. 
4. Equation ( |4.25[ ) can be easily solved by consider the intersections of two functions 
/i(/-0 = ^jr- and ^(a 4 ) = . Note that f2 is concave in the range \i G [l,oo), 

while f\ is a straight line. Thus ( [4.25|) can have at most two real solutions in the 
allowed range. The result is as follows. Below curve I in fig. 1 on the a — b plane, 
which is determined by 



Q'(fj) = 0, and Q"(/x) = 



(4.28) 



13 



The other real solution has g < 1 and so is discarded. 
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( |4.25|) has no solutions in the desired range. Between curve I and the straight line 
(curve III in fig. 1) 

there are two solutions^ 1 < \i 2 < ^3- It can be checked from ( |4.27| ) that /x 2 is a 
maximum of Q with Q"([i2) < while fi 3 is a local minimum with Q"(p3) > 0. At 
curve I, \i2 and fj, 3 merge and move into the complex plane below curve I. As one 
approaches Q129|) from below (c — > 1 + ), [12 — > 1+. Above line c = 1, fi 2 moves outside 



the fi > 1 region to becomes (|4.19|) and Q4-25| ) has only one solution (i 3 . 



5. Below curve III, pi(p 2 ) > |- Above curve III, \x 2 becomes ( |4.19| ) with pi(p 2 ) < \- 
Thus \i2 undergoes a Gross- Witten type phase transition in the large limit. We 
will show in section 5.2 that the transition is precisely the third order Gross- Witten 
transition. 

6. There is an additional curve (curve II in fig. 1), in the a — b plane, determined by 
equation 

Q> 3 ) = 0, and Q( M3 ) = 0, (4.30) 
where the two minima /ii and fis become of equal height 

= V(m,gi) = V(n 3 , g 3 ) < V(fi 2 , g 2 ) • 

Below curve II, one has 

= V(pi, gi) < V(fi 3 , g 3 ) < V(fi 2 , g 2 ) , 

and above curve II 



^3, #3) < = Vinugi) < V(p2,92) ■ 

To summarize, the structure of critical points for ( f4.7| ) in the a — b plane is plotted in fig. 1. 
Below curve I ( |4.28|) , V has a unique minimum ( |4.16|) . Between curve I and line a = 1 
(curve H in fig. 1), there are three critical points pi < P2 < A*3 with /x l5 /j> 3 minima, while 
\i2 a saddle point in the \i — g plane with one negative eigenvalue. At curve I, [12 and \x 3 
merge together. At curve II, \i\ and ^3 exchange dominance and the system has a first 

1 We do not give their explicit expressions in terms of (a,b), since they are complicated and 
not illuminating. We will specify their qualitative behavior below. 
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order phase transition. At curve III, the saddle point ^2 undergoes a Gross- Witten phase 
transition. At the vertical line H, \i\ and \i2 merge together. To the right of line H, hi 
becomes tachyonic, fi2 disappears, and ^3 > 1 remains a minimum. Note that curve II 
and curve III always lie between curve I and line a = 1, but they can be above or below 
each other. Close to a = 1, b = 0, line III lies below curve II and then intersects with and 
rises above it. 

4-3. Critical points for general models 

We now consider the critical point structure of the general effective action ( |4.1| ) and 
( |4.14f) . As is clear from the derivation above, the overall phase structure of Q4.10|) presented 
in fig. 2 to a large extent only depends on the convexity of the function 4 £ and its 
derivative. Our discussion above for (|4.10| ) goes through for a convex S(fx) in (|4.15| ) 
provided that S'(fi) is also convex. For example, for \i < 1, with J-(ij) given by (|4.11|) , 
Q(n) has just one critical point given by the minimum of <S(/i). For \i > 1, the critical 
points of ( |4.15[ ) satisfy the equation (which is a generalization of (|4.25|) ) 

S-W = <i±^ (4.31) 

Since the right hand side of ( |4.31|) is concave and S'(fi) is convex, Q can have at most 
two critical points in the range \x G (1, 00). The pattern and the evolution of the critical 
points with the parameters of the theory also precisely resemble that of ( |4 . 1 (Jj ) including a 
Gross- Witten phase transition for ^2 at (jl = 1. Note that convexity of S'(fj,) implies that 
S'(x) is concave. We thus conclude that the structure of critical points of (|4.1| ) is universal 
if S(x) is convex and S'(x) is concave. In appendix A, we give a more detailed description 
of the phase structure of ( [4.1] ) using the Hartree-Fock approximation. 

4-4- Thermal History 

As discussed in section 3, we would like to use the matrix model ( |3.5| ) as a phenomeno- 
logical model to study weakly coupled string theory in AdS§ x £5 at finite temperature. 
The parameters a, b are functions of the 't Hooft coupling A and the temperature T. In the 
last subsection, we analyzed its large iV critical point structure. In this subsection we will 
show that with some weak assumptions about the A and T dependence of a, 6, the model 
captures all the essential features of the bulk story, which we reviewed in sec. 2. 
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We first identify the critical points of V in the large N limit with the saddle points 
of the Euclidean gravity. \x\ = a has p\ = in large N limit, i.e. the winding in the 
Euclidean time direction is a good quantum number. We thus identify it with the thermal 
AdS background. ^3, which is a minimum, can then be identified with the Euclidean 
big black hole phase (BBH). /i 2 , which has a unique negative eigenvalue can be identified 
with the Euclidean small black hole phase (SBH) in AdS. That a small black hole in AdS 



has a unique negative eigenvalue was pointed out in |31]]. Moreover, one can show from 
the effective potential ( ]126|) and (|4.20| ) that n% has a negative specific heat while (13 has 
a positive specific heat, without knowing their explicit dependence of (a, b) on T. The 
derivation is given in Appendix B. This matches well with the thermodynamic properties 
of the small and big black holes. 

For fixed 't Hooft coupling A, as one varies the temperature T, (a(A, T), b(X, T)) trace 
a curve in the a — b plane, which we will denote by C\. Any curve C\ in the a — b plane which 
starts below curve I in fig. 1 at low enough temperature and ends up to the right of the 
vertical line a = 1 at sufficiently high temperature (assuming it intersects curves I, II, III 
and H only once) reproduces qualitatively the Hawking-Page pictureEl. The thermal 
history following such a C\ can be described as follows. At sufficiently low temperature, 
we start with some point below curve I, where the theory has a unique critical point fj,i, 
corresponding to thermal AdS. As T increases to a temperature To, C\ will intersect the 
curve I (|4.28|) , where new critical points \i<i (SBH) and (BBH) come into existence. At 
Ti > Tq, it intersects with curve II, at which \i\ and \xj, change dominance. This is the 
Hawking-Page transition. Above Ti, 113 (BBH phase) dominates and thermal AdS becomes 
only metastable. As T increases further C\ will eventually hit a = 1 (line H) from the left, 
where the large N expansion (perturbative string expansion) around fii (thermal AdS) 
breaks down. This temperature should be identified with the Hagedorn temperature Tjj 
in AdS string theory. 

There are a few other important features of our matrix model not visible in the bulk 
supergravity analysis: 

1. We found that there exists a line III where the SBH phase undergoes a Gross- Witten 
transition from a gapped phase to a gapless one in the eigenvalue distribution of 
U. We will denote T c the temperature where C\ intersects with line III. Since line 
III lies between line I and line a = 1, we should have T < T c < Tjj. T c can 
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We have drawn such a hypothetic curve in fig. 1. 
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be lower or higher than the Hawking-Page temperature T\ depending on where C\ 
hits line III. This phase transition for SBH, while not visible in supergravity, has a 
natural interpretation in string theory. We would like to identify it with the so-called 



Horowitz-Polchinski correspondence point |?3| for SBH, i.e. the point at which the 
horizon size of SBH is comparable to the string sea In 1 12 311 it was argued that 
as one adiabatically decreases the string coupling, a black hole makes a transition 
to a state of highly excited strings with the same quantum numbers (such as mass, 
charge, angular momentum etc). Here we fix the string coupling (i.e. N), but raise 
the temperature adiabatically. The horizon radius of a small black hole decreases and 



eventually reaches the string scale. Following [23] we would like to argue that beyond 



T c , it is more appropriate to view the critical point (|4.19|) as describing a set of highly 
excited string states. That the phase transition is third order suggests that the energy, 
entropy and specific heat of a SBH vary continuously across the correspondence point 
as it becomes a highly excited string states, but derivatives of the specific heat jump 
in the large iV limit. In the limit of large A, i.e. R 3> Z s , we expect that T c should 
be close to Tjj and much greater than the Hawking-Page temperature T\. That T c 
is below Tjj appears to be consistent with the physical picture of the micro canonical 
ensemble (see e.g. [|32],[nj). For notational simplicity, below we will continue to refer 
to (|4.19|) as the SBH phase, keeping in mind that it should really correspond to a 
highly excited string state. 
2. Our matrix model indicates that at the Hagedorn temperature Tjj , the critical points 
associated with thermal AdS and the SBH merge together. This appears natural since 
very close to the Hagedorn temperature, thermal AdS will be dominated by a few long 
string states. One expects that the distinction between thermal AdS and the highly 
excited string phase which the SBH becomes above T c disappears at the Hagedorn 
temperature. This is again consistent with the physical picture of the micro canonical 



ensemble of [p2| , |10| . 

Above Tjj, i.e. when the C\ curve goes to the right of line a = 1, thermal AdS becomes 
tachyonic and the critical point corresponding to SBH disappears. BBH remains the 
only stable phase. The physical interpretation of a tachyonic thermal AdS is not 
completely clear to us. In the past, it has been argued that string theory above the 
Hagedorn temperature can be interpreted again as some kind of long string phase 



16 Note that the SBH should be considered as a ten-dimensional black hole. 
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which can be analyzed by analytic continuation^! (see e.g. fl551). Our result does not 
contradict this point of view. But we should note that this tachyonic thermal AdS 
phase cannot be reached from the micro canonical ensemble. Whether it plays any role 
in the canonical ensemble is not clear to us. 

We now make some comments on the possible dependence of curve C\ (i.e. curve 
(o(A, T), b(X, T)) with A fixed) on the 't Hooft coupling A. At A = 0, i.e. free theory, 
b(T) = |||10|. Cq moves along the a-axis from a = 0atT = 0toa— > oo at T - > oo. It 
crosses all lines in fig. 1 at a single point (a = 1, b = 0). This was the case analyzed in [HJ. 
At weak coupling it has been shown in P^j that for pure gauge theory at weak coupling, 
b(X,T) = 0(A 2 ) >0anda(A,T) = a(0,T) +0(X 2 ). If the result also holds for M = 4 SYM 
theory, then, Ca<i corresponds to a curve slightly rising above the horizontal a— axis. In 
this case T c < T\. In the supergravity limit A ^> 1, T , T\ ~ while T c , Th ~ f-, i-e. there 
is a big hierarchy between these temperature scales. One can in principle determine part 
of C\ for A — > oo by equating the free energy in the gauge theory with the corresponding 
free energy in supergravity. To leading order in large N we equate the Euclidean actions 
of the SBH and BBH to the corresponding actions of \i2 and \iz in the gauge theory. More 
explicitly, 



^ QM = 4 (4-32) 



where /_ and 1+ are the Euclidean actions (given in (|2.3|) ) of the SBH and BBH respectively 
and Q is given by ( f4.26|) . Since /_ and 1+ are functions of the single variable t = T/Tq, we 
can use ( |4.32|) to determine a(t) and b{t) in the A — * oo limit. Note that this comparison 
is only valid between curve I and curve III, since above curve III, \i2 undergoes a large N 
phase transition which is not visible in supergravity. 

To summarize, our phenomenological (a, b) model reproduces all the important fea- 
tures of string theory in AdS$ x S$ at finite temperature. In fact we got more. We found 
a description of the Horowitz-Polchinski correspondence point in terms of a Gross- Witten 
transition and a non-perturbative picture at and beyond the Hagedorn temperature for 
thermal AdS. 



After finding the critical points one can then use Q4.2|) to find the large N expansion 
of Z(a, b) around them. These expansions should correspond to perturbative string expan- 
sions around the corresponding bulk geometries. In the rest of the paper we will examine 



17 This argument was made in flat space. When the radius of AdS is much larger than the 
string scale, we expect the behavior of thermal AdS to be similar to that of flat space. 
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in detail various regions of the parameter space where these expansions break down and 
study the physics there. Such regions include: 

1. At Tq, where the BBH and SBH saddles merge together. Perturbative string expan- 
sions around BBH and SBH are not valid. 

2. At the Hawking-Page temperature Ti, where there is a first order phase transition 
between thermal AdS and BBH. Although the large N expansion around each critical 
point does not break down, the large N expansion of the full partition function requires 
a special treatment. 

3. At T c , where the Gross- Witten transition for SBH takes place. 

4. At Tjj, the Hagedorn temperature of thermal AdS. 

Note that in 1 and 4 above, the breakdown in the large N expansion happens for the 
metastable phases. In the terminology of the first order phase transitions, To and Th are 
the spinodal temperatures for the BBH and thermal AdS respectively, beyond which the 
metastable phase become unstable. 

4-5. Full partition function and smoothening of Hawking-Page transition 

An immediate consequence of considering the theory at finite N is that the sharp 
Hawking-Page transition is smoothed out to a region of width of order N~ 2 . 
In the infinite N limit, the partition function of the system is 

riogK! + 0(l/A 2 ) T<T X 
\ogZ=l (4.33) 
{-N 2 Q(fi 3 )+ logics + 0(1/N 2 ) T > Ti 

where K\ and K 3 are Gaussian factors computed from the integral ( [4.7| ). Recall that 
Q(Hi) = 0. Qi/J^) equals to zero at T\ and become negative (positive) above (below) T\. 
The transition is first order with a nonzero latent heat given by 



E = N 2dQM 



d(3 

The expectation value of the Polyakov loop also jumps at T\ 



+ 0(1) (4.34) 

Ti 



0(1/A 2 ) T<Ti 
f>((T) = { a (4.35) 

i^ + V 1 -^) T>Tl 
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At finite iV we need to include contributions from both geometries. The full partition 
function of the system between temperature To and Th can then be written in terms of 
the following asymptotic expansion 



Z^ e -^Q(^)A 1 + e- N2 ^A 3 
where A x , A 3 are asymptotic series around fi± and p 3 respectively 

too N 
l + J2N- 2n c n (a,b) 
n=l 
/ 00 

A 3 = K 3 ll + Y,N- 2n d n (a,b) 



(4.36) 



(4.37) 



Note that there is no contribution from ^2, which is a maximum of Q(p)- Below the 
Hawking-Page temperature T x , the p x saddle dominates and p 3 is only metastable. The 
contribution of the second term in fl4.36|) is exponentially small compared with the first. 
Their roles reverse above T x . 

The sharp transition at T x in smoothened out at finite N into a finite region T — T x ~ 
0(N~ 2 ). We now examine this cross over region in some detail. At T y Tl = e C 1, we can 
expand Q(p 3 ) as 



where 



Q(fi 3 (a(T),b(T)), a(T), b(T)) = -vt + 0(e 2 ) 

da 



v = p 2 p 1 + p 4 qx > 



Pi = T x 



dT 



db_ 
dT 



(4.38) 
(4.39) 



and p 2 = p\{ii 3 {T x )). Note the latent heat fl4~34D is related to v by E = N 2 vT x . 
To focus on the transition region we consider the following limit 

N -> 00, e ->• 0, t = eN 2 = finite 

In this limit we find that 

logZ = log {K x (Ti) + K 3 (T x )e vt ) + 0(N~ 2 ) 



(4.40) 



( 4.40|) smoothly interpolates between the first and second line of ( ^4.33| ) as t varies from 
—00 to +00. The expectation value of the Polyakov loop is given by 



p\{t) 



1 dloffZ 



P 



K 3 {T x )e 



vt 



+ 0(1/N 2 ) 



N 2 da r K X (T X ) +K 3 (T x )e» t 
p 2 smoothly interpolates between the first and the second line of (|4.35 ). 
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(4.41) 



5. SBH and tunnelling 

As pointed out in || when T thermal AdS and BBH can tunnel into each other 

with the Euclidean SBH as the instanton bounce [f?3|| . Our effective potential Q(fi) (see 
fig. 2) gives a concrete realization of the physical picture. Through thermal fluctuations, 
BBH or thermal AdS can jump to the top of the barrier to become a SBH. Since it 
has negative specific heat, the small black hole then can either become thermal AdS by 
emitting thermal Hawking radiation or become a big black hole by absorbing radiation. In 
the Euclidean description, this corresponds to rolling down the hill from the two sides of 
the effective potential. In thermal equilibrium, of course the total probability to go from 
thermal AdS to BBH or vice versa should be the same. 



5.1. Tunnelling between thermal AdS and BBH 

In this subsection we will calculate the tunnelling rates between thermal AdS and 
BBH using our effective potential description. For definiteness, we will restrict to the 
temperature range T c > TTq (i.e for (a, b) lying between line I and line III). In this range 
the details of smoothening of singular behavior of the effective potential ( f4.10| ) at \i = 1 
by non-perturbative effects will not be relevant and we will use ( |4.10|) to perform the 
analysis. The discussion for the range Tjj > T > T c is similar, and will not be repeated. 
The only difference is that since the SBH is precisely located at fx = 1 for T# > T > T c , 
it is therefore more convenient to use the two-dimensional effective potential ( |4.8|) . The 
perturbative expansion around SBH breaks down near T c , and will be discused in next 
subsection. 

The tunnelling rates between thermal AdS and BBH can be readily computed using 



the effective potential following the standard procedure p5| , |3q , p7| . To be definite, let us 
first consider the tunnelling rate for thermal AdS over the barrier. In computing ( f4.2| ), 
instead of using the \i contour going from — oo to +oo along the real axis, we consider a 
contour C\ along the real axis from — oo to \xi and then deform the contour at /i2 along 
a steepest descent contour to the complex \i plane (see Fig. 3a). The partition function 
obtained using contour C\ is given by 

Z x « e~ N2Q ^K x (1 + 0(N~ 2 )) + -e~ N2 ^K 2 (l + 0(N~ 2 )) (5.1) 

2 
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In contrast, the flat space has a non-perturbative instability at any finite temperature [34]. 
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where K\ and K<i arise from the Gaussian factor in the saddle point approximation. The 
resulting free energy (let us call it F\) has an imaginary part given by 

I mFl „ ^ e -^v 2 (Q(^)-Q(m))^a (l + (N- 2 )) (5.2) 
The tunnelling rate is then obtained from ( |5.2| ) by |57| 

ri « ^ImFi = ^ e -^ 2 W(M 2 )-Q(m))^ (l + 0(iV- 2 )) (5.3) 

7T 27T ill 

where ujq is the frequency for the unstable mode around the SBH background. Similarly 
the tunnelling rate from BBH to thermal AdS can be obtained by computing ([4.2| ) along 
a contour C2 (see Fig. 3b) to be 



r 2 « ^ e -^ 2 (Q(^)-Q(M3))^ h + (N- 2 )) (5.4) 

27T i\3 




(a) (b) 

Fig. 3: (a) plots the contour C\ which is used to compute ( p.2[ ) and (b) plots the 
contour C2 used to compute (f 



Note that the imaginary part of the partition functions with contours C±, C2 can also 
be obtained by a Borel resummation over the perturbative expansion around the BBH and 
thermal AdS respectively. For example, the asymptotic expansion around the BBH ^3 has 
the form 

e -N*Q{n 3 ) Ks L + \£ N -^ dn {a, b)j (5.5) 
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Due to the presence of a maximum at \x 2 , the coefficients d n at large n are given by 

K 2 Y{n) 



The asymptotic expansion ( |5.5| ) is clearly divergent and not Borel summable. In fact the 
Borel transform dr £ of -J^ is singular at z = Q(/j> 2 ) ~ Qi^), thus preventing 
Borel summability if the difference between the action of SBH and BBH Q(n<2) — Qil^s) 
is positive. This is a familiar situation in instanton physics, where we can interpret the 
singularity in the Borel transform in the positive axis as real instantons. One can integrate 
the Borel transform just above the singularity on the real positive axis, producing an 
imaginary part. Using this procedure, we find that ( |5.5| ) becomes^ 

e-N'Q^Ks (1 + 0(N~ 2 )) - - e ~ N2 ^K 2 (l + 0(N~ 2 )) (5.7) 

which is precisely what we get by computing the partition function following contour C 2 in 
Fig. 3b and leads to the tunnelling rate ( |5.4|) . Similar discussions apply to the perturbative 
expansion around thermal AdS. Note that the instanton effect we have obtained associated 
with the SBH is of order e 9 ? for g s = i and not of order e~~^ typical of D-instantons. 
In some sense, the effect we are describing could be interpreted as a collective state of N 
D-instantons. 

Of course the total partition function is real and therefore cannot contain any imag- 
inary part. Physically this is equivalent to a balance rule between the probability to 
nucleate a SBH from thermal AdS and the probability of decay of the metastable BBH. 
The two imaginary parts that we obtain in contours C\ and C 2 should cancel in the com- 
plete partition function that is real. This is in fact the case once we realize that the sum 
of contours C\ and C 2 gives the real contour which defined the total partition function Z. 

To conclude this subsection we will briefly comment on the asymptotic expansion for 
T < Tq. Below line I the two saddles \i 2 and fis move into the complex plane, in other words 
the real instanton that we have associated with the SBH above line I becomes complex. 
This in particular means that the corresponding singularity in the Borel transform is not 
any more on the real positive axis and Borel summability is potentially restored^. From 
the physical point of view this means that the thermal AdS saddle is stable below line I. 



19 One keeps only finite number of terms in the first term below. Also there is a sign ambiguity 
here. We take the sign to be the same as that of contour C2. 

20 This phenomena is well known in quantum field theory. For instance for the X(p 4 theory 
we have complex instanton solutions that produce a priori harmless singularities in the Borel 
transform on the negative real axis. 
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5.2. A Gross-Witten transition for small black hole 

In this subsection we examine in detail the behavior of SBH as it crosses line III 
(c = 1) in the a — b plane. It undergoes a third order phase transition in the large N limit 
and the large N expansion around SBH breaks down there. We will define a double scaling 
limit to smooth out the transition. 

When c = 1 + e with < e « 1, we find from (f4.23|) - (|425|) that the critical point 
corresponding to the SBH is given by 

e 2 e 1 e 

^ = 1 + - + ..., £ 2 = i + _ + ..., /9l (« 72 ) = _ + _ + ... (5.8) 

When c = 1 — e, one finds from (|4 .19 ) that 



02 = 1, 2 = l-£ + ..., pi (g 2 ) = __! + ... (5.9) 

Thus as c crosses 1, gi crosses 1 and the eigenvalue distribution of 02 crosses p\ = \. One 
can also check that the second derivatives around the SBH varies smoothly across c = 1. 

As we commented after equation ( fll2| ), (|4.8|) has a third order discontinuity at g = 1 
which is smoothened out at finite N. More precisely, let 

g = l-N~*y (5.10) 

Then V should be replaced by0 

1 2 i 00 

" 71 = 

where the F n (y)'s are smooth functions of y. In particular, Fq describes the doubling 
scaling limit of (|4.4|) with the following asymptotic expansion 

£ - ± log(-y) + ■ ■ ■ -y > 1 

/'n(//) = { , x (5.12) 

Note that ( ^.11|) smoothly interpolates between the second and the third line of (4.8). 
21 See e.g. sec. 3.1 of [19]. F n below correspond to there. 
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Equations (|5.8|) - (|5.11| ) suggest that to study the physics of the SBH near c = 1, we 
can consider the following scaling 

a(T) = a + ai eq, b(T) = b + hqe, 2 ( 1 ~ a ") = 1 ( 5 . 13 ) 

bo 

H = l + e 2 x, g = l-ey, e = iV~§ (5.14) 

where ao = a(T c ), ai = T c a'(T c ) (similar for b) and eq = T 7 r Tc . Thus we have 

2(1 — a) ai bi 1 , n . . 

c = = 1 - cieq, ci = h — = 7-( 2a i + Oi) 

6 1 - a Oo oo 

Note ci > according to our assumptions on the thermal history of the model and q < 
as we approach line III form below. 

Plugging ( |5.13|) and (|5.14|) into (|5.11j ), we find that 



» 2 v= 1A ^ L + - 2 *(y-'f)-My) + o(e) (5.15) 

Note that in the last expression, x simply plays the role of a Lagrange multiplier. The 



integral over x,y in (|4.7|) is not well defined since we are integrating around the neighbor- 
hood of a saddle point in the x — y plane. If we rotate the integration contour of x to be 
along the imaginary axis, the x integral will result in a delta function for y and we find 
that the partition function around the SBH is given by 

/7T JV 2 (l-a) 2 . „ /cig\ 
_ e i^+M—) (1 + 0(e)) (5.16) 

The prefator i can be understood as due to the tachyonic mode of the SBH. We find 
that around T c , one can define a double scaling limit where the SBH is described by 



Fq. It was argued in |38[] that F (t) describes the full partition function of the type OB 
theory in d = dimension, i.e. pure 2-d supergravity. The parameter t is proportional 
to the cosmological constant \x in the super-Liouville interaction. We are then led to the 
conclusion that in a double scaling limit around T c (as we argued earlier in the Horowitz- 
Polchinski correspondence point) the SBH appears to be described by type OB theory in 
zero dimension. 
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6. Catastrophes and the break down of perturbative string expansions 

In this section we examine the critical behavior of the BBH at Tq where its saddle 
merges with that of the SBH and the critical behavior of the thermal AdS at the Hagedorn 
temperature Tjj- We find in both places that the breakdown of the large N expansion can 
be understood in terms of the simplest type of catastrophes allowed by the symmetry The 
divergences at the perturbative level can be smoothened out at finite N using the standard 
techniques of catastrophe theory 

6.1. Nucleation of black holes 

The large iV expansion around the big black hole saddle breaks down near line I, 
where it coalesces with the unstable small black hole saddle. We will show that the critical 
behavior there is given by the fold catastrophe. One can define a double scaling limit in 
which the partition function for this sector is given by an Airy function. 

From QjjgD , curve I can be parameterized by 



l-2w 

a[w) 



(l-w) 2 (l + w) 



2w 

b ( w ) = T\ \2ft i_ >>3 ' W e [°> X ] 

(1 — W) Z {1 + w) 6 



(6.1) 



Suppose that at temperature Tq, the curve (a(T),b(T)) intersects the curve I ( |6.1| ) at a 
point labelled by wo, i.e. a (To) = o(wo), b(T Q ) = b(wo). At the intersection point, ( |4.25| ) 
has a double root given by wo, which is an inflection point of Q{w). We will consider the 
behavior of Q(w) ( |4 . 2 6| ) near w = wq and T = Tq. Let 

a(T) = a(T ) + a x e, b(T) = b(T ) + he, w = w + y (6.2) 

with e = T y T ° and ai = T a'(T ), b\ = T b'(T ). We will consider the regime |e| 1 and 



\y\ ~ V l e l ^ 1- Plugging (|6.2| ) into ( |4.26| ) and expanding in e and y we find tha 

Q = C (e, w Q ) - f (~y 3 + qey) + 0(e 2 , y\ y 2 e) (6.3) 



22 The expansion below breaks down at wo = 0, where it can be checked that Co and various 
higher order terms become singular. At Wo = 0, we have ao = 1, bo = and the physics goes over 



to that of |J. 
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with 

/= -2 aW ' 9= iHrr°-J (6 - 4) 

where tan6>o is the slope of line I at wq. Note that 

a'(w ) < 0, b'(wo) > 0, tanflo = V/^"! = ~7TT~ — v?> 6> G (tt/2, tt) 

From our assumption of the thermal history of the theory, q is positive. Note that Co(e, Wo) 
is analytic in e. 

For e > 0, from ( |6.3| ) Q has one maximum and one minimum at y — ±^/eq. The two 
extrema merge at q = and move to complex values for e < 0. The values of Q at the 
minimum and the maximum are given by (e > 0) 

Qo = C (e)T^ 1 e^+0(e 2 ) (6.5) 

The second term gives the leading nonanalytic term in e and the specific heat has a critical 
exponent 7 = i. 

It is instructive to compare ( |6.5| ) with the result ( |2.3| ) from supergravity. Expanding 
( EOp in e around (3q for big and small black holes we find 

I = C (e)T2ei + ••• (6.6) 

with Co an analytic function of e. We see exact agreement between ( |6.5|) and ( B.6|) in the 
critical exponent. The critical exponent i is universal, depending only on the fact that a 
maximum and a minimum merge together (fold catastrophe). 

As T — > To, the large N expansion around the big and small black holes break down. 
The physics around them can be captured by a double scaling limit. Since we are interested 
only in the BBH and SBH, we will again consider ( |4.2|) along contour C2 in Fig. 3b. From 
(|T3|) we introduce a new variable 

z = qe(N 2 f)i 

and consider the scaling limit 

e — » 0, iV — >• 00, z = finite 
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In this limit, the partition function becomes 



Z 2 =C I dse-^ s3+zs (6.7) 
Jc 2 

which is given by an Airy function. Cq is a non-universal factor 

where T\ is the 0(1/N 2 ) term in J 7 (fx) (|4.11| ) (not given explicitly there). It is easy to 
check that our choice of contour gives the Airy function 

Z 2 = 2niC e-^ 1 Ai(ze-^ 1 ) (6.8) 

( |B~8|) smoothes out the divergences in perturbative expansions. Also note that the argu- 
ment of the Airy functions precisely sits on the Stokes line. This is a consequence of the 
fact that for ^2,^3 real, ImQ((i 2 ) = ImQ^) = 0. Although ( |6~8| ) is complex, the full 
partition function should be real when including the contribution from the contour C\ of 
fig. 3a.i 

6.2. Hagedorn behavior for thermal AdS 

We will now examine the merger of thermal AdS and the SBH (or more precisely long 
string phase). For this purpose, let us first look at the free energy near the thermal AdS 
background for a < 1, which can be found by expanding the partition function around the 
saddle (|4.16| ). Around ( f4.16| ), the partition function can be approximated as 



N f 1 N^ ( „ ,2 1 

Z X = ^= dfxe- — ^-^- — 



- a) 



N 



00 r°° N r 
/ dxe~~^ 

n J —00 



2Vnb 



(l_ a )n+i (6.9) 



^T(n + i) ( 2Vb ^ n 



1 -°^o ^ \N(l-a) 

There are also nonperturbative corrections which are omitted here. Thus the free energy 
around the thermal AdS background can be written as 

2b 106 2 2965 3 

logs, = - log(l - «) + ^ + + 3(1 - a)>N> + ' ' ' (fU0) 



23 We would like to thank G. Festuccia and A. Scardicchio for extensive discussions regarding 
this point. 
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Note that in contrast to the free theory case analyzed in [19], the free energy now receives 
perturbative contributions to all orders. The free energy diverges as a — > 1. The leading 
order divergence, arising from genus one contribution, is 

logZi = - log(T - T H ) + const (6.11) 



since as a — > 1, 1 — aocT — Tjj- Q6.ll ) is precisely the Hagedorn divergence for string 



theory in a spacetime with all directions compactified (recall that AdS behaves like a box). 
By a Laplace transform of ( |6.11| ) one finds that the density of states is given by 

tt(E) w const e pnE (l + 0{l/E 2 )) (6.12) 

Also note that as T — > Tjj , the free energy around SBH is given by 

logZi oc -N 2 (T-T H ) 2 + ■■■ (6.13) 



The perturbative expansion (|6.10|) breaks down at 1 — a ~ N 1 . To explore the 
physics near this point, we let0 

a = l-N~ 1 q, n = l-N~ 1 x, g = 2N~^y (6.14) 

and consider a double scaling limit with N — > oo with q,x,y finite. We find that 

Q = N- 2 P + 0(N~ 3 ) 

with 

P= (a: ~ g)2 + X y 2 . (6.15) 
4o 

Note that P has two critical points for q > 



(1) *i = g, yi = 0, P-=(» °) 

(2) x 2 = 0, y 2 = A /^, ' ^ 



26' \ /2a 

b 







2 Below for notational simplicity we will assume b does not change as a = 1 is crossed. To 
incorporate the change in b is straightforward. 
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corresponding to (fii,gi) and (1^2,92)- They merge at q = and for q < only the first 
solution remains. The relevant part of the partition function ( f4.7|) then becomes 



at r°° r°° 
Z x = -= dx ydye- p{x ' v) 

V Tib J -00 Jo 

= 2N J ydye- qy2+by4 



(6.16) 



Introducing a (O-dimensional) complex scalar (p with y = \<p\, the second line of ( |6.16| ) 
can also be written as 

Z x = J cU/xfy* exp [-m 2 ((3)(j)*<P + b(<P<P*) 2 ] (6.17) 

with m 2 (/3) = q oc T — Tjj. It is tempting to identify <p with the so-called thermal scalar 
(a winding tachyon) in string theoryll. Indeed ( |6.17|) coincides with the effective action 



one expects for a thermal scalar near Hagedorn temperature ||39|| . It is clear from the 



second line of (|6.16|) or ( |6.17|) that the merger of thermal AdS and SBH is described by a 



cusp catastrophe with q = (i.e. a = 1) corresponding to the cusp point. Note that the 
appearance of effective action fl6.17|) is forced on us by the U(l) symmetry of the complex 



scalar field <p, which corresponds to the symmetry of the boundary effective action in 
the large iV limit. The cusp catastrophe is the simplest possibility consistent with this 
symmetry. The only nontrivial dynamical input in ( |6.17| ) is that b > 0, following from the 
existence of a first order Hawking-Page transition. 

The integral in ( |6.16| ) is not bounded as is expected since we are focusing in the 
neighborhood of the effective potential containing only thermal AdS and SBH. To define 
the integral in ( |6.16| ) we will choose an integration contour analogous to C\ of Fig. 3a in 
section 5. We take the contour in the y-plane to go from y = to the maximum 1/2 = 
along the real axis and then go straight up to the complex value at y2- The integration 
along the real axis will give us an error function, which smoothes out the divergences of 
the perturbative expansion at q = 0. 



25 Note that we see a zero dimensional scalar since AdS may be considered as a box. 

32 



7. Conclusions and discussions 



In this paper we introduced a class of phenomenological models to understand string 
theory in AdS 5 x S5 in the canonical ensemble. Our models reproduces all the known 
qualitative features of the theory. They also have some interesting predictions includ- 
ing the existence of a third order phase transition for SBH, which we identify with the 
Horowitz-Polchinski point. We studied the simplest model ( |3.5|) in great detail. We found 
the Hagedorn behavior of thermal AdS at Th and the critical behavior of nucleation of 
Euclidean SBH and BBH at T are governed respectively by cusp (A 3 ) and fold (A 2 ) 
catastrophe. It is clear these features persist for all models in the class due to universality 
of the catastrophe. We believe this gives strong indication that they capture qualitative 
behaviors of a weakly coupled string theory in a large AdS spacetime. Since for a large 
radius AdS, the SBH resembles a ten dimensional Schwarzschild black hole in flat space- 
time, and the Hagedorn behavior for strings in AdS resembles that of flat space, we expect 
that the behaviors we observe here may yield clues to answers for similar questions in flat 
spacetime. 

There are many other questions which can be explored along the lines of our investi- 
gation. For example, it would be nice to have a better understanding of our proposal that 
the Horowitz-Polchinski point for a small black hole should correspond to a Gross- Witten 
transition in the boundary theory. It would be interesting to understand whether the 
process involves changing the spacetime topology. Also, since the saddles corresponding 
to thermal AdS and SBH merge at Th, we expect that in the worldsheet sigma model 
of thermal AdS, turning on the marginal operator corresponding to the thermal scalar at 
Th, the theory can be deformed into a SBH background. It would also be interesting to 
understand from the worldsheet point of view what happens when the tachyon in thermal 
AdS above the Hagedorn temperature condenses^. The discussion of ||0| might be useful 
for this purpose. 

We believe that the phenomenological approach developed here can have many other 
applications. For example, it would be interesting to develop an effective potential ap- 
proach for the tunnelling discussed in . It would also be interesting to see whether one 
can use our methods to address the problem of black hole information loss. 



26 It appears clear from the effective potential picture that the theory will flow to the BBH 
background. 
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Appendix A. Large N phase structure and thermal history for general matrix 
model flOD 

For completeness, in this section we discuss the phase structure and thermal history 
of matrix model (|4.1|) using the Hartree-Fock approximation. The Hartree-Fock treatment 
of double trace operators was earlier discussed in |[42|| . The Hartree-Fock approximation 



gives equations of motion which can be solved to find the critical points of the theory and 
the value of the action evaluated at the critical points. However from this method one 
cannot find the off-shell effective potential (essential for our purposes). 

In the infinite N limit, it is convenient to introduce the density of eigenvalues 

1 - 

P(«) = ? E^- ff i)' ~7V<0<7V (A.l) 



with 



( |4.1| ) can be written as 



N 

i=l 



±TrU n = p n = I (Wp(O)i' 



Z = J[Dp]e- h TV W (A.2) 

where V[p] has the form 

V[p} = -± J ded<f>p(0)p(<t>) Plog\2 S m -^) -S(\ Pl \ 2 ) (A.3) 
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Since the potential is symmetric, we can take p\ to be real. The equations of motion 
following from (|A.3| ) can be written as 



/ 



Q X 

dfj>p{(j>) cot — — - = Ksm6 (A.4) 



with 

k = 2S'(x)pi , x = pi . 

The solutions to the above equation are well known [ pO] , |2"T| , |2"2|] , leading to the self-consistent 
equations for p\ (using ( }4.6| ) with g replaced by k above) 



(A.5) 



pi = S'(x)p 1 , < pi < ^ 

S ' W = 4p 1 (l 1 - Pl )' U?^ 1 
Note that the first line of ( |A.5| ) implies that 

Pi = (A.6) 

or 

S'(x) = l, x = pje [0,1/4] (A.7) 



We can slightly rewrite the second equation of (|A.5|) as 

S'(x) = f(x), f(x)= 1 xe [1/4,1] (A.8) 

4 v /x(l - ^x) 

Note are also functions of 't Hooft coupling A and temperature T. We will show 
that given the following assumptions about S(x; A, T), (|A.2| ) has exactly the same large N 
phase structure as that of the (a, b) model analyzed in the main text and thus that of the 
AdS supergravity: 

1. S(x) is convex, i.e. S'(x) is monotonically increasing; 

2. S'(x) is concave; 

3. For sufficiently low temperature, S'(x) lies below f(x) defined in (|A.8|) in x G [1/4, 1]. 
For sufficiently high temperature S'(0) > 1; 

4. 5"(l/4) is a monotonically increasing function of T. 



We note that for (a, b) model (|3.5| ), conditions 1 and 2 are automatically satisfied. Con- 
dition 3 corresponds to our assumption that the C\ curve starts below line I of fig. 1 at 
sufficiently low temperature and ends to the right of line H at high enough temperature. 
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Condition 4 makes sure that C\ intersects all lines in fig. 1 only once as T is varied. Note 
that conditions 4 may be further relaxed. 

We first note that f(x) is a monotonically increasing and convex function which takes 
values 

/(1/4) = 1, /(l) = +oo (A.9) 

At a sufficiently low T, condition 3 implies (|A.8| ) has no solution. Since S' (1/4) < 1 and 5" 
is monotonically increasing, equation ( |A.7| ) also does not have a solution. The only phase 
of the system at this temperature is thus given by ( |A.6|) , i.e. 



xi = 



This is the thermal AdS. It is a minimum of V for S'(0) < 1. As we increase the tem- 
perature, the curve S'(x) will start intersecting^ with f(x). The temperature at which 
they become tangent is To, where the critical points corresponding to SBH and BBH start 
appearing. Immediately above T , ( |A.8| ) will have two solutions | < X2 < X3 < 1 and from 



condition 2 it can only have two solutions. Since V is bounded from below, X3 should be 
minimum (BBH), while X2 a maximum (SBH). At this temperature ( |A.7|) again does not 
have any solution since S' (1/4) < 1. At a temperature T c > Tq, when S' (1/4) = 1, X2 = \ 
is both a solution of flOp and ( gg ). When T > T c , S"(l/4) > 1, flA~8|) only has one 
solution X3. X2 moves to the region x < 1/4 and becomes a solution to ( |A.7| ). Convexity of 
S(x) implies that the solution to ( |A.7| ) is unique. At a temperature Th, when S'(0) = 1, 
X2 and x\ coincide. Above Tjj, %2 no longer exists and x\ becomes tachyonic due to that 
S'(0) > 1. Since at Th, x\ and X2 coincide, we have V(x2) = V(xi) = in the large N. 
Since V(x2) > ^(^3), ^(^3) must be smaller than zero at Tjj. We thus conclude that 
there must be a first order Hawking-Page transition at some temperature Tq < T\ < Th- 
To summarize, the general model satisfying the four assumptions above have exactly 
the same phase structure as that of the (a, b) model including the Gross- Witten phase 
transition for SBH and the merger of SBH and thermal AdS at Th- 



27 This is guaranteed following our assumptions. 
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Appendix B. Specific heat of small and big black holes 

In this appendix, we show that the phases corresponding to \ii and fis have negative 
and positive specific heat respectively. 

We first show that (|4.19|) has a negative specific heat. From ( |4.20| ), 

,d 2 V 



2 n2 ' 



-N z p 



dp 2 



N 2 p 2 (da 1 -a db 



2b \df3 b d/3 



(B.l) 



< 



We now look at the specific heat of a solution to equation ( |4.25|) . Evaluated at a 
critical point /x c (/3), the action is given by Q(f3, fi c ((3)). We first note an identity 

2 



d 2 



dp 



d 2 Q 
dfj, 2 



d 2 Q 



dp 2 



(B.2) 



In deriving the above equation we have used the equation of motion ^ | Mc = and 

d 2 Q 



dfJ-c 

d/3 



d 2 Q 



(B.3) 



We note from (|4.26| ) that 



d 2 Q 



dp 2 
d 2 Q 



dpdn 



1 f da n c — a db 
2b \d[3 + b dp 



1 f da /i c — a db 
~2b \d/3 + ~~b 



dp 



(B.4) 



Thus the specific heat for \i c can be written as 

c v {nM) = -N 2 p 2 -^Q{P,n c {p)) 
N 2 P 2 /da + fic- a db 



Ab 2 \dp 



dp 



1 



d 2 Q 
dfi 2 



(B.5) 
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For /j, c = /j<2, §^|„„ < 0, and we find 



H2 



c v ((i 2 ) < . 



For n c = it follows from p.27| ) that 

o<Q"M = -k ' 



l 



Plugging ( p.6| ) into ( P-5| ) we find: 



2b V 4v/^T> 3 /2 2b 



c v (fi 3 ) > . 



(B.6) 
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